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Abstract 

We establish correspondances between factorisations of finite abelian 
groups ( direct factors, unitary factors, non isomorphic subgroup classes ) 
and factorisations of integer matrices. We then study counting functions 
associated to these factorisations and find average orders. 
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1 Introduction 

In this paper we study the correspondance between different factorisations of fi- 
nite abelian groups and that of integer matrices. In particular, we are interested 
in the matter of enumeration. This is in line with the questions on direct and 
unitary factors of abelian groups first raised by Cohen |Coh60j in 1960 and sub- 
groups of abehan groups studied since long in different contexts (see, eg, Butler 
|But87j . Goldfeld-Lubotzky-Pyber |C;LP()4j or Bhowmik-Ramare |BR98j i. 
We consider the decomposition of a finite abelian group G as a direct sum 

G = Z/niZ © ■ ■ ■ © Z/nin2'L © ■ ■ ■ © Z/nin2 ■ ■ ■ rir-Z 

where rij are divisors of the order n = n\rC^^ ■ ■ - rir and r the rank of the group. 
In terms of integer matrices this group can be associated to a canonical represen- 
tative of a double coset 

M{G) = diag[ni, nin2, ■ ■ ■ , nin2 ■ ■ - nr]- 
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The matrix M{G) is said to be in Smith Normal Form (SNF) or sometimes an 
elementary divisor and its determinant equals the order of the group G. It is 
well known that the number of non-isomorphic abelian groups of order upto x 
is asymptotic to Lx |ES34j where L is a constant. If the rank is bounded by r, 
this corresponds to the number of r x r SNF matrices of determinant upto x and 
then the constant depends on r |Bho93j . 



In 1960 E. Cohen |Coh60j defined the factorization of finite abelian groups as a 
decomposition in formal direct factors. He introduced the zeta function associated 
to the direct factorization, G = H (B K, and obtained that the number of direct 
factors, on an average, is log x. In Section 1 we interpret direct factorisations of 
G in terms of decompositions into SNF factors of the 'conjugate' matrix of G. 
The idea of conjugate matrices comes from considering the type of a group, in 
the sense of Macdonald |Mac79j . Thus we can treat the truncated zeta function 

C^(.)C^(2.)c^(3.)--.c^M 

as a counting function of certain direct factors of abelian groups. The case of 
unitary factorisations is treated similarly. 

Later, |BR98J it was proved that a subgroup of G corresponds to a divisor of 
M{G) which is a canonical representative of a right coset, for example a matrix in 
Hermite Normal Form (HNF). Using the zeta function associated to subgroups it 
was found that on an average there were many more subgroups than direct factors. 
The number of subgroups of an abelian group is essentially x^^^ which 
is heavily dependent on the rank. 

In Section 2 we consider classes of subgroups under isomorphism. In terms of 
matrices we are interested in divisors of M{G) which are in SNF. We now inves- 
tigate whether on an average the number of non-isomorphic subgroups depends 
on the rank of the group. 

To do so we introduce the zeta function 



n)n 

n=0 

where Ar{n) denotes the number of subgroup classes (under isomorphism) of 
abelian groups of order n rank atmost r and use this function to get the mean 
value 

X 

^^ylr(ra) ~ CrX log X 

n=l 

for positive integers r. This shows that on an average the number of non- 
isomorphic subgroups is about the same as the number of direct factors whereas 
the number of all subgroups is much larger as the rank increases. The same 
average order holds true when we consider subgroups of arbitrary rank. 
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While Zr{s) in general have no global functional equation we will prove that 
for all positive integers r it has a local functional equation 

Zr{s)= n Fr{p-'), 
p prime 

with 

This follows from a recursion formula that we develop for the functions Fr{x). We 
include some values of the generating function by using this recursion formula, 
the first few being 

Zi{s) = C (s), Z2{s) = — — and ^3(5) = G{s), 

where G{s) is a Dirichlet series absolutely convergent for Re(s) > 1/5. 

For the limit case Z{s) = lim^^oo Zr{s) we show that while Z{s) does have a 
meromorphic continuation to Re(s) > 0, it does not have a meromorphic contin- 
uation beyond the imaginary axis. We tend to believe that the same is true for 
Zr{s) when r > 3. 

Acknowledgement The second author thanks Radoslav Dimitrig for useful dis- 
cussions. 

2 Direct Factors. 

In this section we discuss the relationship between direct factors of abelian groups 
and factorisation of matrices. We recall a definition. 

Definition 1. An abelian group G has a (formal) direct factorisation into H and 
noted asG = HQ)KiiG = H + K,H and K being subgroups of G whose 
only intersection as subgroups is trivial. 

Formal direct factors of abelian groups were studied in the context of arith- 
metic functions of abelian groups and Cohen |Coh60j defined direct Dirichlet 
convolution of abelian groups in terms of this decomposition, i.e. 

(Xi*X2)(G)= Yl XiiH)x2{K), 

G=H®K 

where Xi ^^re arithmetical functions of abelian groups, in an attempt to generalise 
the factorisation of integers. 
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With no loss of generalisation, we can consider p-groups, i.e. 

Gp ^ Z/p^'Z © Z/p^'^Z © • • • © Z/p^'-Z, 
where the order of Gp is p«i+«2+- +ar^ prime number p and < oi < 02 < 

• ■ ■ < ttr- 

The invariants of G is given by the partition A = (ai, 02, • • • , a^), which, in 
combinatorics, is called the type of G. We call the group whose type A' is the 
conjugate of A, the conjugate of G and denote it by G'. We now express the 
direct factorisation above as a matrix factorisation. 

Definition 2. The SNF factorisation of a SNF matrix 5" is a decomposition into 
SNF matrices whose product equals S. 

We now establish the following bijection 

Proposition 1. If M{G') = M{H')M{K') is a SNF factorisation, then G = 
H ® K is a direct factorisation and conversely. 

Proof. Let G be a group of type 

A=(^, /l + /2,---, /l + /2 + --- + /m> 
gi times 92 times Qm times 



Then 



A' = ( Qm , 9m + Qm-l; ■ ■ ■ , 9m + Qm-l H Vgi) . 



Writing gi = hi + ki, sum of non-negative integers, we get types of groups H' and 
K' such that M{G') = M{H')M{K') is a SNF factorisation, i.e. 

M{H') = dia9 \p^^, . ^. . . . . . 

fm times fm-l times /i times 

etc. whereas 

H ^ ®hi^/pf'Z Z/pf'+f^Z ■■■®hm Z/pf'+f^+-+f"^Z 

and 

K ^ ©fc^Z/Z^Z ®k, Z//i+^^Z ■■■®km Z//i+^2+-+-^'"Z, 

where ©^ before a summand means that it is repeated t times. 
Thus the type of H is 

/l + /2,..., /l + /2+--- + /rn ) 
hi times ft.2 times hm times 
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and that of K is 



/l + /2,---, /l + /2 + 

ki times k2 times 



times 



and the relation G = H Q) K is satisfied. 

The converse is proved by taking G 
steps. 



H (B K as above and retracing the 

□ 



Example 1. Let G = Z/p'^Z © Z/p'^Z © Z/p^Z. Here A = (2,2,4), thus the 
conjugate A' = (1, 1, 3, 3). The following non-trivial SNF factorisations of M{G'), 
i.e. 



M{G') 



/I 0^ /p 

10 

p 

\0 pj 



0\ 

p 

p2 

\0 p2y 



/I 

1 

p2 

\0 






py 



fp 0\ 

p 

p 

\0 p/ 



give H = Z/p^Z,K = Z/p'^Z © Z/p^Z and H = Z/p^Z © Z/p'^Z,K = Z/p^Z. 
Obviously, changing the order of the matrices would interchange H and K. 



The generating function for the number of direct factors of a finite abelian 
group is given by 

as was first done by Cohen |Coh60j and subsequently studied by many, for ex- 
ample Kratzel |Kra88j . Menzer |Men95j . Wu |WuOOj or Knopfmacher |Kno85j . 
A truncated form of this function, 

occurs as the Dirichlet series associated to the SNF factorisation of an r xr integer 
matrix, see, for example |Bho93j . which we can now interpret in the context of 
abelian groups: 

Theorem 1. The function ('^{s)('^(2s)('^{3s) ■■■ ('^{rs) generates the direct fac- 
tors of finite abelian groups whose conjugates have rank atmost r. 

This implies the following result on average orders: 

Corollary 1. The number of direct factors of non-isomorphic abelian groups of 
order atmost x whose conjugates have rank atmost r is asymptotic to A^xloga;, 
where is a constant that depends on r. 
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3 Unitary Factors. 



Unitary factors of finite abelian groups were also studied by Cohen as a general- 
isation of the corresponding idea on integers, where one encounters the function 

t{n) = E 1- 

d\n, gcd{d,n/d)=l 

Definition 3. An abelian group G has unitary factors H and K if G = H (B K 
is a direct factorisation H and K have no further common direct factor. 

As above, we see that in terms of matrices this corresponds to block SNF 
factorisations, i.e. for M{G') as above, M{H') is a unitary factor if hi is either 
or Qi for each i. Thus the number of unitary factors of G is 2"*. 

Example 2. With G as in Example 1, the only non-trivial unitary factorisation 
H = Z/p^Z ® Z/p^Z, K = Z/p^Z corresponds to the SNF factorisation 

Ph \ ^ //2 \ fph 
p3j J I p2j J I pi, 



M{G') -- 

where I2 is the 2x2 identity matrix 



The level function Ur{n) that counts the number of unitary factors of abelian 
groups of order n whose conjugates have rank atmost r is given by 



where uj{n) gives the number of distinct prime factors of n. We obtain 
Theorem 2. The associated zeta function 



n=l 

can be written as 

Uris) - 



C(2s)C(4s)C(6s)---C(2rs) ' 

This implies the following result on average orders: 

Corollary 2. The number of unitary factors of finite abelian groups of order 
atmost X whose conjugates have rank atmost r is asymptotic to Krxlogx, where 
Kr is a constant that depends on r. 

In the infinite dimensional case, the study of the average order first occurs in 
literature at the same time as direct factors and continues to be refined. (See, 
eg. Calderon |(]ain3j or Zhai [ZhaOnp . 
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4 Subgroup Classes 



It is easy to see that the number of classes of subgroups of G corresponds to the 
number of divisors of M{G) which are themselves in SNF. Thus we are interested 
in the function 

E i = Ei. 

Mi|M(G), Ml SNF H<G 

where for all isomorphic subgroups H, we count only once. For our purpose it is 
enough to consider a p-group 

Gp ^ z/|j"iz e • • • © z/p"i+-+"'-z 

whose type is the partition (ai, ai + 02, . . . , ai + • • • + a^), > 0. We wish to 
count all subgroups H whose type is "less than" the type of G, i.e. where the 

type of H is {hi, 61 + 62, • • • , &i H VK) with XlLi ^» ^ XlLi 1 < J < ^■ 

Thus we study the function 

ai a-i+a.'z—hi ai+a2H \-ar—bi—h2 h, — i 

/^(ai,a2,...,ar) = ^ XI " ' XI 

61=0 62=0 hr=0 

Collecting the exponents roi + • — \- ar — we write 

00 

Fr{x) ^Y.°'r{k)x\ (2) 

fe=0 

where 

ar{k) = ^ /r(ai,a2, . . . ,a,.). (3) 

raiH |-ar=fe 

Now let Cj = oi + • • • + Oj and (ij = 61 + • • ■ + hj and let = (ci, C2, ■ • • , c^) 
and = ((ii, (i2, ■ ■ ■ , ^r) be partitions of k and m <k respectively. Then ar{k) 
denotes the number of pairs of partitions {(Afe,/im)}, such that di < q. 

We first obtain a recursion formula for the ar{k) . To do so directly seems 
difficult, and it is convenient to introduce an additional parameter j. 

Definition 4. Let ar{k,j) be defined as ar{k), but where the condition < di < 
■ ■ • < dr is replaced hj —j < do < di < ■ ■ ■ < dr, i.e. 

Cl C2 Cr 

M*,.)= E E E E-- E 1- w 

Cl-\ \-Cr=k do=—j di=do d2=dl dr=dr—l 

0<Cl<-<Cr 
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Proposition 2. One has the recursion formula 

""<*'^' = |o. k > 1. 

iii) arik, -j) = 0, j > 1. 

[k/r] 

(Hi) arik,j) = arik,j -l) + '^ar-iik-mr,m+j), (r > l,j > 0). 

m=0 

(iv) arik, 0) = ar(fc)- 

Proof. Equations (i), (ii) and {iv) follow from Definitional 

For {Hi), the condition —j < do in Definition 4 can be divided into two cases, 
either — (j — I) < do from which we get the contribution ar{k,j — 1), in {iii) ; or 
do = —j in which case we let Ci = m and denote 

c- = Cj+i — m, and c?^ = di+i — m. 

It is clear that (0, c'^^, C2, ■ ■ ■ < c'^_i) is a partition of — mr., while d'^ < 
0, — m — i < d'^ < d\ < ■ ■ ■ < d'^_^, and c?^- < c'j. We thus get a 
contribution 

ar^i{k — mr, m + j) 

for each < m < [k/r], and this concludes the proof of {iii). □ 
We introduce the notation a{k) = ak{k). 

Since the number of partitions of k into k non-zero parts is equal to the 
number of partitions of k into r parts, r > k, we get, 

ar{k) = a{k). (r > k) (5) 

Thus, 

lim ar{k) = a{k) 

We now obtain some simple bounds in terms of the classical partition function. 
Lemma 1. We have that 

q{k) < a{k) < kq{k)^, 
where q{n) denotes the total number of partitions of the n. 
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Proof. The lower bound follows from the definition. Now let = (ci, C2, ■ ■ ■ , Ck) 
and \m = {di, d2, - ■ ■ , d^) be partitions of k and m < k respectively. If di ^ 0, 
can be completed to (1, ■ ■ ■ , 1, d\, d^, - ■ ■ , cifc), a partition of k. There are atmost 
g(fc) of such Am- 

If di = 0, 7^ 0, we can complete Am to get a partition oik — c\. In this way, 
for each Am we have atmost 

q{k) + - Ci) H V q{k - Ck) < kq{k) 

possibilities. Since there are q{k) ways to write A^, the upper bound follows. □ 

We will now turn our attention to the generating functions 

Definition 5. Let ar{k), ar{k,j), and a{k) be as above. We define the generating 
functions to be 

oo 

Fr{x) = y^^ar{k)x'', 

k=0 

oo 

Fr{.x,y) = ^ ar{k,i)x^y\ 

k,j=0 

and 

oo 

F{x) = 

fe=0 

These definitions are consistent, since 



Frix) = Fr{x,0), (6) 

the only term that survives with y = is when j = 0. 

Lemma 2. With F and F^ defined as above one has that F{x) and Fr{x) are 
analytic functions in the unit disc with integer power series coefficients such that 
F{0) = Fr{0) = 1. Furthermore the function F satisfies the inequality 

Fix) > (0 < X < 1) 

Proof. The power series coefficients of F^ and F are integers since they are count- 
ing functions and by Proposition 2 (i) and (iv) and eq. we have that ar(0) = 1 
and a(0) = 1, which implies Fr{0) = F{0) = 1. By the well known generating 
function for the classical partition function 

oo ^ 
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and the lower bound in Lemma ^ 

q{k)<a{k), 

this gives us the lower bound in Lemma |21 Equation (|7j) also implies that the 
generating function of the partition function is analytic in the unit disc, and hence 
the classical partition function q{n) is of subexponential order. This implies that 
k{q{k)Y is of subexponential order and by the upper bound in Lemma El so is 
a{n), and also ar{n) since < ar{n) < a{n). This proves that F and are 
analytic in the unit disc. □ 

We will now see how the recursion formula for the coefficients ar{k,j) in 
Proposition 121 yield a recursion formula for its generating function Fr{x,y). 

Proposition 3. One has the following recursion formula for the function Fr{x, y) : 

1 



(0 Fo{x,y) 
(ii) Fr{x,y) 



Fr-i{x,x'')x'' -Fr-i{x,y)y 



{x-" -y)(l-y) 
Proof. From Proposition |21 (ui) it is clear that 

oo [k/r] 

Fr{x,y){l - y) = ar~i(fc - mr,m + j)x''y\ 

k,j=0 m=0 

and by replacing k and j suitably, we get 

oo j 

Fr{x,y)-Fr{x,y)y= ^ ^ x^+™V-™«r-i(A:, j), 

k,j=0 m=0 

k,j=o V y / 

_ Fr-i{x,x'')x'' - Fr-i{x,y)y 
x'^ — y 

Proposition |21 (i) and Definition 5 give 

oo 

i=o 

1 



□ 
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We first notice that 
Lemma 3. The function F satisfies the identity. 

Proof. This follows from eq. (jH} and Proposition IHl (m). □ 

We now calculate the first few values of Fr{x). 
Lemma 4. One has that 

(t) Fo{x) = 1, 

and 

(m) F2{x) = (1 - x)-\l - x'^)-^{l - x^). 

Proof. Putting ?/ = in Proposition El («) proves (i). From Lemma El we get that 

Fi(x) = Fo(x,x) 



which proves (ii). By Proposition 121 (m) we get that 

p. X Foix,x)x - Fo{x,y)y 

Fi{x,y) = r- ^ 

{x-y){l-y) 

xy — 1 

By Lemma El we get 

x^ — 1 



(8) 



F^ix) = F,{x,x^) 



(X-I)2(x2-1)3' 

which proves (Hi). □ 
We further prove: 

Proposition 4. The functions Fr{x, y) and -F^(x) satisfy the functional equations 

{^) F^ix,y)=x-^^^^'y'y-'FJ-,-), 

\x yj 

Frix) = X-'^'+'^/'Fr (^i^ . 
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Proof. We will use induction to prove (i). Let r = 0. Proposition El («) gives us 
that 

-2 



1 



y 

X y 



1 1 



and hence (i) is true for r = 0. We now assume that (i) is true for r = k. By 
using Proposition (jH)) we now obtain 

Using the functional equation (i) for r = A; it equals 

^ -(fc+l)(fc+2)/2^2 -^A:(a^ \x ^' ^ — \y ^)|/ ^ 

= x-^'^'^^'^'^''y-'F,^,{x-\y-\ 

which finishes the proof for (i). We will now prove (u). By Lemma El we have 
that 

Fr{x) = Fr-l{x,x'), 

The functional equation with y = x'' gives 



X 



\x x^ J 
= x-'^'+'^/^F^ ( -] . 



X 



□ 



5 The zeta function of subgroup classes 

5.1 The zeta function of subgroup classes of bounded rank 

We will now introduce the zeta function. We use the notation Ar{n) for the 
multiplicative function that counts the number of isomorphic subgroup classes of 
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abelian groups of rank less or equal to r and order n. If n is a prime power then 

Ar{p^) = ar{k). 
The corresponding zeta function is given by 
Definition 6. Let 



oo 



n=l p prime 



We then get 
Lemma 5. One has that 



m=l 



where Gk,r{s) is a Dirichlet series without real zeroes and absolutely convergent 
for Re(s) > l/{k + 1). Furthermore the I3r{m) 's are integers. 

Proof. By Lemma 2 the result follows from Dahlquist |Dah52j with 

Mm) = J2^,[^)^Br{d), (9) 

d\m 

where B,. comes from the expression 



logF,,(a;) = Y,Br\ 



m)x^. 



m=l 



Thus 

fe 

p prime m=l 
k 

m=l 

is without real zeroes and absolutely convergent for Re(s) > l/(fc + 1). □ 

With the notation /3(m) = we get from eq. (0) and the definition of 

the I3r{m) eq. ^ that 

(3{m) = (3r{m). (r > m) (10) 
By using Definition 6, Lemma eq (fTUI) and Lemma 4 we obtain 
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Theorem 3. One has that 

Zi{s) = C{s), 

and 

Zkis) = eis)e{2s)C'{3s)Gs,k{s), (k > 3) 

where 6*3,^(5) is a Dirichlet series without real zeroes and absolutely convergent 
for Re(s) > 1/4. 

This implies the following result on average orders: 

Corollary 3. The number of subgroup classes of abelian groups of rank less than 
or equal to r and order atmost x is asymptotic to CrX log x where Cr is a constant 
depending on r. 

5.2 The zeta function of subgroup classes of arbitrary 
rank 

We now let A{n) = linir^oo ^r(^)5 Z{s) = \imr-,oQ Zr{s) and proceed as in the 
last section. We obtain with Lemma 2 and Dahlquist's theorem that 

Lemma 6. One has that 

k 

Z{s) = n amsf'^"'^Gk{s), 

m=l 

where Gk{s) is a Dirichlet series without real zeroes and absolutely convergent for 
Re(s) > l/(fc + 1), and the [3{m) 's are integers defined by eq. (fTUj) . 

which implies 

Theorem 4. Let A{n) denote the number of subgroup classes of abelian groups 
of order n, and let 



Z{s) = Y,A{ 



n]n 



n=l 



denote the corresponding zeta function. Then 

Z{s) = eisK\2s)CH3s)G,is) 

where 6*3(5) is a Dirichlet series without real zeroes and absolutely convergent for 
Re(s) > 1/4. Furthermore Z{s) is a meromorphic function for Re(s) > and 
the imaginary axis is the natural boundary for Z{s). 
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Proof. That Re(s) = is the natural boundary follows from Lemma 2 which 
shows that F can not be written as a finite product Y['j=ii'^~^'')"^^ ■ ("^i ^ ^) 
and from Dahlquist's theorem the function cannot be meromorphically continued 
beyond the line Re(s) = 0. □ 

As before we get an estimate: 

Corollary 4. The number of subgroup classes of abelian groups of order almost 
X is asymptotic to Cxlogx where C is a constant. 

Remark 1. The constant C in Corollary 4 can be calculated by 

oo 

C = JJ C(^)^^'"^ = 13.1854452968422695. 

m=2 

We notice that j3{m) > for m = 1, . . . , 12 and for A; = 13 Lemma El allows 
us to see (with the Mathematica program in the Appendix) 

Remark 2. One has that 

xa7s)C\8s)C\9s)C\l0s)e\l2s)iai3s))-'Gn{s), 

where Gi3(s) is a Dirichlet series absolutely convergent for Re(s) > 1/14 In 
particular this means that there are no poles for Z{s) for Re(s) > 1/13 except 
for s = 1/k, k integer. Under the Riemann hypothesis, we can use Lemma 6 
with = 25 to see that there are no poles for Z{s) for Re(s) > 1/26 except for 
s = 1/k, k integer. 

For more calculations of /5(m) see the appendix where we use a Mathematica 
program to calculate its value for the first 100 values of m. 

5.3 A question on the zeta function of subgroup classes 
of bounded rank 

It seems likely that r = 1 and 2 are the only cases where Zr{s) has a meromorphic 
continuation to the entire complex plane. For higher rank we expect that 

Zr{s) is a meromorphic function for Re(s) > with the imaginary axis as its 

natural boundary. 

That Zr{s) has a meromorphic continuation to Re(s) > of course follows from 
Lemma 5. To prove that Re(s) = is the natural boundary we would have to 
show that Fr{x) cannot be written as a finite product 

k 

Fr{x) = JJ(1 - x^r^. {mj e Z) 
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Thus it is sufficient to show that not all roots of the numerator of Fr{x) he on 
the unit circle which is easily verifiable by numerical calculation for small values 
of r. As an example, let r = 3. We can use the explicit form of ^3(2:) as given in 
the appendix, and see that the equation 1 + 2x + 2x^ + + = Q has a root 
X = -0.621744 + 0.440597i with an absolute value 0.762031. 
Such calculations are not obvious for a general r. 

6 Appendix 
6.1 Tables 



Table of (3r{m) 



r\m 
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-6 


6 


1 


9 


4 
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4 


-4 
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2 


9 


-10 


-2 


-2 


3 


25 


-25 


11 


5 


2 


3 


3 


4 


2 


-1 


-6 


-1 


3 


4 


2 


3 


-12 


-22 


15 


26 


6 


2 


3 


3 


4 


2 


6 


-7 


-4 





-2 


1 


13 


10 
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-4 


-9 


7 


2 


3 


3 


4 


2 


6 


1 
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-5 


12 


-1 


8 


17 


-9 


8 


2 
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3 


4 


2 


6 


1 


4 


-4 


-8 


-8 


6 


-2 


-3 


18 


20 


9 


2 


3 


3 


4 


2 


6 


1 


4 


6 


-9 


-11 


3 


-8 


-4 


7 


21 


10 


2 


3 


3 


4 


2 


6 


1 


4 


6 


2 
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2 


6 


1 


4 


6 


2 





12 


-15 


-16 


-3 


3 


13 


2 


3 


3 


4 


2 


6 


1 


4 


6 


2 





12 


-1 


-17 


-6 





14 


2 


3 


3 


4 


2 


6 


1 


4 


6 


2 





12 


-1 


-2 


-7 


-3 


15 


2 


3 


3 


4 


2 


6 


1 


4 


6 


2 





12 


-1 


-2 


9 


-4 


16 


2 


3 


3 


4 


2 


6 


1 


4 


6 


2 





12 


-1 


-2 


9 


13 



For the case of P{m) we have the first 100 values: /5(1), . . . , /9(100) = 2, 3, 3, 4, 
2, 6, 1, 4, 6, 2, 0, 12, -1, -2, 9, 13, 0, -16, 6, 44, -25, -6, 16, -19, 52, -21, 52, 
-103, -140, 505, -203, -381, 286, 88, 185, -751, 564, 1015, -2304, 1007, 1876, 
-3432, 1177, 3665, -1582, -6119, -2558, 21792, -7745, -34936, 40625, 4248, 
-34948, 19176, -361, 30668, -104511, 81530, 196372, -457425, 148194, 497952, 
-459549, -261973, 163195, 1015378, -808365, -1895457, 3467251, -483924, 
-4303223, 4879125, 535751, -4991204, -634706, 7525125, 7551313, -31891601, 
10748210, 53812155, -65485063, -16092470, 74326725, -12532241, -51617936, 
-53724288, 219808510, -91391802, -403364243, 663219049, -34111622, 
-908127917, 645523852, 643374980, -402710571, -1752056675, 1385661209, 
3632651805, -6352203822, -317398867 
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6.2 Mathematica Programs 

We can now implement the recursion formulae in Mathematica. 
Clear [F] ; 

F[0,x_,y_] := F[0,x,y]=l/(l-y)-2; 
F[r_,x_,y_] := F[r,x,y] = 

Factor[(x^r F [r-1 .x.x^r] -F [r-1 ,x,y] y)/( (x"r-y) (1-y) )] ; 

Now using 
F_r(x)=F[r,x,0] 

gives us for r = 0, 1, 2, 3 that 

Fo{x) = 1, 

FAx) 



1 



F,{x) 



-l + xf' 
1 + ,x + 



3 ' 



and 



Fsix) 



-1 + x) {1 + x) 



{l + 2x + _^ ^3 _^ ^^^^ (1 + X + 2x^ + 2x^ + x"^) 
{-l + xf{l + xf{l + x + x^)'' ' 



The program above is slow if we just want to calculate the lower order coefficients 
ar{n) in the power series of Fr{x), so the program we used to calculate the Table 
of I3r{m) is the following (in this case with m = 16) : 

Clear [a] ; 

a[r_,n_,-l] :=a[r,n,-l]= 0; 

a[0,n_,j_]:= a[0,n,j] = If [n==0, j+1 ,0] ; 

a[r_,n_,j_] :=a[r,n, j]=a[r,n, j-l]+Suin[a[r-l,n-m r,m+j] , 

{m.O, Floor [n/r] }] ; 
a[r_,n_] :=a[r,n]=a[r,n,0] ; 

m=16; Table [aa=Series [Log [Sum[a[r , j] *x" j , {j ,0,in}] ] , {x,0 ,m}] [ [3] ] ; 
Table [ b=Divisors [n] ; 

Sum[MoebiusMu[b[[j]]]/b[[j]] aa[[n/b[[j]]]] , {j , 1 , Length [b] }] , 
{n.l.m}] ,{r,l,m}] 

A similar program is used to calculate the /5(m) for m = 1, . . . 100. 
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